Dominance of excitation continuum in the longitudinal spectrum of weakly coupled 

Heisenberg 5 = 1/2 chains. 



O 

o 

(N 

> 

o 
O 
in 

(N 



A. Zheludev, 1 '^ K. Kakurai, 2 T. Masuda, 3 K. Uchinokura, 3 and K. Nakajima 4 

1 Solid State Division, Oak Ridge national Laboratory, Oak Ridge, TN 37831-6393, USA. 
2 Advanced Science Research Center, Japan Atomic Energy Research Institute, Tokai, Ibaraki 319-1195, Japan. 
Department of Advanced Materials Science, The University of Tokyo, Tokyo 113-8656, Japan. 
* Institute of Solid State Physics, The University of Tokyo, Tokai, Ibaraki 319-1106, Japan. 

(Dated: February 1, 2008) 

A low-field spin flop transition in the quasi one-dimensional antiferromagnet BaCu2Si2C>7 is ex- 
ploited to study the polarization dependence of low-energy magnetic excitations. The measured 
longitudinal spectrum is best described as single broad continuum, with no sharp "longitudinal 
mode", in apparent contradiction with the commonly used chain-MF/RPA theories. The observed 
behavior is also quite different than that previously seen in the related KCUF3 material, presumably 
due to a large difference in the relative strength of inter-chain interactions. The results highlight 
the limitations of the chain-MF/RPA approach. 
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The physics of weakly-interacting spin chains is of 
great fundamental interest, as it represents the crossover 
from quantum to semi-classical spin dynamics. Exci- 
tations in conventional three-dimensional magnets are 
single-particle spin wave states. These can be described 
as precessions of magnetic moments around their equi- 
librium orientation, and are necessarily polarized per- 
pendicular to the order parameter. In contrast, the 
excitation spectrum of onc-dimensioanl (ID) quantum- 
disordered antiferromagnets (AFs) is isotropic. In par- 
ticular, in the 5=1/2 Heisenberg model, the spec- 
trum is a polarization-independent multi-spinon contin- 
uum CD, |^, |3|, |], Js). Arbitrary weak inter-chain interac- 
tions in quasi-lD 5 = 1/2 AFs restore long-range order 
(LRO) §, 0, |, g, break rotational invariance and rein- 
state the transverse spin waves at low frequencies [ll] . 
The continuum persists in the high-frequency regime. 

A counter-intuitive novel feature of weakly-interacting 
Heisenberg 5 = 1/2 chains is the longitudinal mode (LM), 
a long-lived "spin wave" polarized parallel to the direc- 
tion of ordered moment Jl^, [ll[] . This new type of exci- 
tation was recently discovered in the KCUF3 compound 
[fl2| , [l3|] , where it is seen as a single symmetric broad peak 
in the longitudinal spectrum. |l2j Such behavior is not 
fully understood within the commonly used chain-Mean- 
Field H and Random Phase Approximation (MF/RPA) 
theories [[[o| [ll] that predict a sharp (long-lived) LM 
and a separate broad longitudinal continuum (LC). To 
further probe the limitations of the MF/RPA approach, 
and learn more about the exotic LM, we performed a 
polarization-sensitive inelastic neutron scattering study 
of BaCu2Si2C>7. The advantage of this material is a con- 
siderably smaller ratio of inter-chain and in-chain interac- 
tions, as compared to KCUF3. This brings us closer to the 
idealized weak-coupling limit, where the MF/RPA theory 
can be expected to work best. Surprisingly, the measured 
longitudinal spectrum substantially deviates from theo- 
retical predictions. Contrary to expectation, the LM in 
BaCu2Si207 is altogether ill-defined, and the low-energy 



longitudinal spectrum is best described as a single broad 
continuum feature. Thus, for longitudinal polarization, 
there is no analogue of the separation of single-particle 
and continuum states, that was previously seen in the 
transverse spectrum jl4[ [TH) . 

The magnetism of BaCu2Si2 0y is by now very well 
characterized using bulk methods 0, |l7|], neutron 
diffraction fill 



18 and inelastic neutron scattering |L 



|15| , [L6J, [19]. The material is orthorhombic (space group 
Pnma, a = 6.862 A, b = 13.178 A, c = 6.897 A) 
with slightly zigzag AF 5 = 1/2 chains of Cu 2+ ions 
running along the c axis. The in-chain exchange con- 
stant is J = 24.1 meV. Interactions between the chains 
are much weaker, and the characteristic badwidth of 
spin wave dispersion perpendicular to the chain direc- 
tion ("mass gap" in the terminology of Ref. |0|) is 
A = 2.51 meV. BaCu2Si207 orders antiferromagnetically 
at Tjy = 9.2 K= 0.033 J/Ur with a zero-T saturation mo- 
ment of only mo = 0.15 /ib- For comparison, for KCUF3, 
J = 37 meV, T N = 39 K= 0.09J//c B , m = 0.54 ^ B 
(Ref. H H; |H] ancl references therein). Due to a weak 
easy-axis anisoltropy in BaCu2Si2C>7, the ordered stag- 
gered magnetization is aligned along the c axis. 

The most direct way to determine the polarization of 
magnetic excitations is by using polarized neutrons. Such 
measurements are typically associated with severe inten- 
sity losses in the polarizing monochromator and analyzer. 
Given the size of currently available single-crystal sam- 
ples, and the expected scattering cross section for longi- 
tudinal excitations in BaCu2Si207, this type of experi- 
ment does not appear feasible. Fortunately though, even 
for an unpolarized neutron beam, the scattering inten- 
sity is polarization dependent. Only spin components 
perpendicular to momentum transfer contribute to the 
cross section. One commonly used strategy is to measure 
the inelastic signal at equivalent wave vectors in differ- 
ent Brillouin zones (BZs). Comparing such data sets can 
provide polarization information, but relies heavily on 
the exact knowledge of the spectrometer resolution func- 
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tion. The latter can be estimated quite accurately, but 
an error of only a few percent can ruin the polarization 
analysis. This was the main problem in previous incon- 
clusive attempt to apply this technique to longitudinal 
excitations in BaCu 2 Si207, as described in Ref. [ p^[ . 

In the present work the problem was overcome by using 
the "reverse" approach. The same inelastic scan (always 
in the same BZ) was performed for different orientations 
of the ordered moment in the system. Polarization in- 
formation was then extracted in a point-by-point data 
analysis that is robust and insensitive to any resolution 
effects. As described in detail below, to change the di- 
rection of ordered moment we exploited one of the re- 
cently discovered held-induced spin-re orientation transi- 
tions in BaCu2Si2C>7 A 5 mm diameter and 40 mm 
long cylindrical single crystal sample was mounted on 
the the thermal 3-axis spectrometer PONTA (5G) at the 
JRR3-M reactor at Japan Atomic Energy Research In- 
stitute at Tokai. The sample cylinder axis is close to the 
(1, 1,0) reciprocal-space direction and was oriented ver- 
tically in the experiment. The (0, 0, 1) and (—0.5, 1.86, 0) 
reciprocal-lattice vectors determined the horizontal scat- 
tering plane. Sample environment was a horizontal-field 
6 T cryomagnet. The field was at all times applied along 
the c axis of the crystal. The data were collected in an 
energy scan close to the ID AF zone-center (0, 0, 1). The 
design of the horizontal-field magnet imposes stringent 
constraints on experimental geometry, allowing only nar- 
row ±15° windows for the incident and scattered beams. 
For this reason, each data set was collected while simulta- 
neously scanning the energy transfer hui and momentum 
transfer perpendicular to the chain axis q± , at a constant 
momentum transfer along the chains q» = 2tt/c, which is 
the ID AF zone-center. Figure |l| shows the trajectory of 
the scan relative to the previously determined spin wave 
dispersion in BaCu2Si2C>7. Also shown is the evolution 
of the instrument resolution ellipsoid, calculated in the 
Cooper-Nathans approxomation J2l], Q . 

The scan was repeated in three independent measure- 
ments. The first measurement (Fig. 0, solid circles) was 
performed in a magnetic field H\ = 1.2 T at T = 2 K. 
Under these conditions, just as at H = 0, the ordered 
moment in BaCu2Si2 07 is parallel to c jl8[ |l9). The 
measured intensity is related to the longitudinal (paral- 
lel to the ordered moment) and transverse (perpendic- 
ular to the ordered moment) dynamic structure factors 
S^{ti,u)) and S x (q.,w) through Ji(q,w) cx S , - L (q,w)(l + 
cos 2 a) + 5" (q, uj) sin 2 a + B(q,w). Here a is the angle 
between the momentum transfer q and the chain axis, 
and S(q, u>) is the background, assumed to be of non- 
magnetic origin. In our particular case a < 15° through 
the entire scan, and it is mostly the transverse spin 
fluctuations that are observed. In the a second mea- 
surement the same scan was repeated at H2 = 2.2 T, 
which exceeds the critical field H c = 2.0 T of a spin-flop 
transition |L7|, |l8). These data are plotted with open 
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FIG. 1: (E — q)-space trajectory of the inelastic scan and 
the corresponding evolution of the FWHM resolution ellip- 
soid (shaded ellipses) . The solid curve shows the dispersion of 
spin waves in BaCu2Si2C>7. The dotted line and light shaded 
area represent the dispersion of the longitudinal mode and 
the domain of the excitation continuum in BaCu2Si2C>7, re- 
spectively, as predicted by the MF/RPA model. 



symbols in Fig. g. In the spin-flop state the ordered 
staggered moment is perpendicular to the field direc- 
tion and therefore lies within the crystallographic (a, b) 
plane. The measured inelastic intensity can be written as 
7a(q,w) cx S ,± (q,w)(l+sin 2 a)+5 ll (q,cj)cos 2 Q:+6(q,w). 
For small a one observes an equal mixture of longitudi- 
nal and transverse spin fluctuations. In the final mea- 
surement, performed at zero field, the sample was re- 
moved from the magnet and sample holder. The same 
scan was repeated to directly measure the extrinsic back- 
ground B, mainly due to spurious scattering in the sam- 
ple environment (triangles in Fig. |J). The intrinsic back- 
ground (from scattering in the sample itself) was previ- 
ously shown to be negligible At each (q, w)-point 
the equations for 1\ and I2 were then solved to extract 
the dynamic structure factors for different polarizations, 
as plotted in Fig. |[ 

Using this method to determine the polarization of spin 
fluctuations is based on several assumptions. The first is 
that the spin structure rotates as a whole when going 
through the spin-flop transition, but is otherwise unaf- 
fected. For BaCu2Si207, such behavior was confirmed 
by previous neutron diffraction studies in magnetic fields 
up to 5 T fjg| ]. In particular, upon crossing if c , the mag- 
nitude of the ordered moment remains almost constant, 
and the change in relative alignment of nearby spins is 
insignificant. The second requirement is that magnetic 
anisotropy and Zeeman effects are negligible on the ex- 
perimentally relevant energy scales, as defined by the 
gap A and the experimental energy resolution (typically 
1.7 meV FWHM). In our measurements this condition 
is also satisfied. Indeed, at H = H2 = 2.2 T the Zee- 
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FIG. 2: Inelastic scans measured in BaCu2Si207 near the 
ID AF zone-center (0, 0,1) at T = 2 K in magnetic fields 
applied along the c-axis. The data were collected below (solid 
circles) and above (open circles) the spin-flop transition at 
H c = 2.0 T. The background (triangles) was measured as 
described in the text. The solid line is a Lorentzian fit to the 
background scan. 

man energy per spin is an order of magnitude smaller: 
gSjj/^H ~ 0.25 meV. Performing differential measure- 
ments at Hi and if 2, with gSfiB^H of only 0.13 meV, 
further minimizes the effect, as shown by arrows in Fig. [ij. 
The magnitude of magnetic anisotropy in BaCu2Si207 is 
also negligible, the associated spin wave gap being only 
about 0.2 meV Jl5| . In summary, the spin flop transition 
only changes the coordinate system in which the terms 
"longitudinal" and "transverse" are defined, but should 
not affect the nature of the excitation spectrum. 

Previous measurements of transverse spin fluctuations 
111 BaCu 2 Si 2 7 @ were found to be in remark- 
ably good quantitative agreement with MF/RPA theory 
JTo| , O ]. In this approach each individual spin chain is 
viewed as subject to an effective staggered exchange field 
Htt, generated by LRO neighboring chains. The stag- 
gered field in turn generates an attractive potential be- 
tween spinons, producing 2-spinon bound states. Two 
components of the spin-triplet bound state are polar- 
ized perpendicular to H^. These are analogues of con- 
ventional spin waves, and have a gap A oc ff^ /3 . The 
third component is the longitudinal mode with a gap 
A|| = V3A. Regardless of polarization, a multi-magnon 
continuum, the remains of the spinon-continuum, has a 
gap of A c = 2 A. Starting from this MF solution the 
RPA is used to calculate the spin wave dispersion per- 
pendicular to the chains. For BaCu2Si207 the MF/RPA 
dynamic structure factor for transverse spin waves is ex- 
plicitly written out in Eqs. 4-7 in Ref. fL5fl . The contin- 
uum part is well approximated by the "truncated Muller 
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FIG. 3: Transverse- (a) and longitudinal-polarized (b) dy- 
namic structure factor measured in BaCu2Si207near the ID 
AF zone-center. The heavy solid line is a fit based on 
MF/RPA theory, as described in the text. The grayed and 
hatched areas represent the single-mode and continuum parts 
of the cross section. In (b) the heavy dotted line is a fit that 
does not include any single-mode contribution, as described 
in the text. The thin dash-dot line is the experimental energy 
resolution. 

ansatz" functional form, as given by Eqs. 10 and 11 in 
the same reference. As an important consistency check, 
this model cross section function, after a numerical con- 
volution with the spectrometer resolution, was fit to the 
transverse spectrum measured in this work. Since all pa- 
rameters of the model were previously determined with 
very good accuracy [ ft5| , |l9[ ] , the only parameter varied in 
the present data analysis was an overall intensity scaling 
factor. An excellent fit was obtained, and is shown in a 
solid line in Fig. ^a. Continuum and single mode parts 
are shown a hatched and greyed areas, respectively. 

We now turn to discussing the main result of this work, 
namely the measured longitudinal spectrum. As shown 
in Fig. ||b, the scattering is a broad asymmetric peak with 
a maximum around 4 meV energy transfer, and an ex- 
tended "tail" on the high-energy side. As a first step, we 
analyzed the data using a model cross section based on 
MF/RPA theory, as for transverse excitations. The cor- 
responding dynamic structure factor for the LM is given 
by Eqs. 12 and 13 in Ref. Jljj]. At the particular wave vec- 
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tor, the MF/RPA predicts the LM at 4.2 meV, quite close 
to the observed intensity maximum. Given the shape of 
the observed peak, in our fits the LM was allowed to 
have a non-zero intrinsic energy width. The correspond- 
ing (5-function in Eq. 12 of Ref. was replaced with 
a Lorenzian profile of FWHM T. In addition, the pref- 
actor 7 that chracterizes the energy-scaled intensity of 
the longitudinal mode relative to that of transverse spin 
waves was treated as an adjustable parameter. Contin- 
uum scattering was assumed to be exactly as for trans- 
verse polarization. This model, when convoluted with 
the spectrometer resolution function, yields a reasonably 
good two-parameter fit to the data with 7 = 2.2(1) and 
T = 2.4(7) meV (heavy solid line in Fig. |b). 

The obtained parameter values clearly indicate the fail- 
ure of the MF/RPA model to describe longitudinal ex- 
citations in BaCu2Si20y. First, the refined intensity co- 
efficient 7 is about 4.5 times larger, than the MF/RPA 
prediction ||rj| . Second, to obtain a good fit, we had to 
assume a substantial intrinsic width for the LM. In con- 
trast, the LM in the MF/RPA model is infinitely sharp. 
For weakly-coupled chains like those in BaCu2Si20y, this 
damping effect has a particularly important consequence: 
the LM becomes altogether ill-defined. Indeed, in this 
case, the continuum is relatively strong, and comparable 
in intensity to single-mode scattering. In the MF/RPA 
the LC has a (pseudo)gap at 2A, and the LM is cen- 
tered at about 1.73A. At the same time, our analysis 
shows that in BaCu2Si2 0y the intrinsic FWHM of the 
LM is itself be about A. Given these numbers, the LM 
and the LC are too close to be separated, and consti- 
tute a single continuum feature. To illustrate this fact, 
in Fig. ||b we also show a fit to the longitudinal data us- 
ing only the truncated Muller ansatz (heavy dotted line) , 
with no single mode contribution. With a continuum gap 
A c = 1.4(1)A (smaller than the expected MF/RPA value 
of 2A), this empirical function describes the data at least 
as well as the cross section that also includes a damped 
LM. The situation is quite different for materials with 
inter-chain interactions of intermediate strength, such as 
KCUF3. Here the LC, which is the high-energy "tail" 
seen in constant-g scans, is relatively weak. As a result, 
even though the LM is broadened, it is seen as a well- 
defined symmetric inelastic peak ||l2| . 

Comparing the spectra measured in BaCu2Si207 and 
KCuF 3 , we conclude that the MF/RPA picture based 
on distinct LM and LC contributions does not become 
more accurate in the weak-coupling limit. The width 
of the LM peak normalized by A in the weaker-coupled 
BaCu2Si207 system is even greater, than in the KCUF3 
compound. This counter-intuitive behavior may be re- 
lated to the critical nature of spin correlations in isolated 
5=1/2 chains. Quantum criticality implies that there 
is no characteristic energy scale in the chains, and there- 
fore any inter-chain coupling is to be considered as strong. 
Strictly speaking, the RPA is well justified only for inter- 



chain interactions small compared to the energy gap in 
individual chains. In our case, the chains are intrinsically 
gapless. The gap A only emerges in the MF theory and 
is necessarily of the same order of magnitude as inter- 
chain coupling itself []l0|]. The problem is therefore only 
marginally suited for the RPA for intrinsic reasons. 

In summary, despite a small ratio of inter-chain and 
in-chain interactions, the low-energy longitudinal spec- 
trum of BaCii2Si207 is dominated by a broad excitation 
continuum, and there appears to be no distinct "longitu- 
dinal mode" . Further theoretical calculations beyond the 
MF/RPA are needed to better understand this behavior. 
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